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We study proximity effects close to a boundary between s and s± superconductors. Frustration, 
caused by interaction of the s-wave gap parameter with the opposite-sign gaps of s± superconductor, 
leads to several anomalous features. In the case of strong frustration a nontrivial time-reversal- 
symmetry breaking (TRSB) state, with nonzero phase angles between all gap parameters, is possible. 
In a more typical state, the s-wave order parameter is aligned with one of the s± gaps. The other 
(anti-aligned) gap induces negative feature in the s-wave density of states, which can serve as 
a fingerprint of s± state. Another consequence of the frustration is an extended region in the 
parameter space in which s-wave superconductivity is suppressed, despite being in contact with 
nominally stronger superconductor. This negative proximity effect is always present for the TRSB 
state, but extends even into the aligned states. We study these effects within a simple microscopic 
model assuming dirty limit in all bands, which allows us to model the system in terms of minimum 
number of the most relevant parameters. The described anomalous features provide a route to 
establishing the possible s± state in the iron-based superconductors. 



I. INTRODUCTION 

The discoverjffi of the iron-based high-temperature su- 
perconductors has brought new theoretical and experi- 
mental challenges to condensed matter physics. Despite 
the undoubted progress in our knowledge (see, for ex- 
ample, Refs. I2H11) some of the most intriguing questions 
still do not have a definite answer. First among them by 
importance is the precise form of the superconducting or- 
der parameter. The most plausible candidate so far is the 
extended s-wave or s± stat^. In this state the gaps on 
the (hole-like) bands at the center and the (electron-like) 
bands at corners of the Brillouin zone (BZ) have opposite 
signs. The physical origin of this state is in the repulsive 
interband interactions, which likely dominate the pair- 
ing in these materials. Although unconventional, such 
state from symmetry point of view is indistinguishable 
from a conventional s-wave state, as they both belong 
to the Aig representation of the lattice rotation group. 
Existence of such state is supported by theoretical cal- 
culations done within the framework of several methods 
- Random Phase Approximation (RPA)^^, Functional 
Renormalization Group (FRG^P^, Fluctuation Exchange 
(FLEX^iS!^ as well as analytic one-loop RG and diagram- 
matic 

calculationpHia. 

In spite of a concentrated effort, the structure of 
the order parameter in these superconductors has not 
yet been unambiguously established from experiment. 
The angle resolved photoemission spectroscopy (ARPES) 
reveals uniform gaps on different bandP^ for several 
compounds close to optimal doping, but it cannot re- 
solve the most crucial issue - the relative sign of the 
gaps in the electron and hole bands. This sign can 
be probed by some phase-sensitive experiments, simi- 
lar to ones performed for the cuprate high-temperature 
superconductor^l^. Even thoug h suggestions for such ex- 
periment have been madd^^'^ they have not been real- 
ized in practice yet.l^^ 



The strongest support in favor of the s±-state comes 
from inelastic neutron scattering experiments'^*^, which 
detect the emergence of a resonant magnetic mode be- 
low the superconducting transition, as expected for such 
sign-changing state. This mode was detected in almost 
all iron based superconductors and its frequency scales 
approximately proportional to transition temperatur^. 
However, straightforward interpretation of the data is 
complicated by the multiband character of the Fermi sur- 
face (FS) and the possible strong role of interactions. 

Indirect probe of the bulk order parameter struc- 
ture is provided by the low-temperature behavior of 
the thermodynamic and transport properties which is 
sensitive to the presence of quasiparticle states at the 
Fermi level. Extensive studies have demonstrated a very 
rich behavior; all compounds fall into three relatively 
well-defined groups. Several clean materials (LiFeAs^il, 
Bai_a:Kj.Fe2As2^^) do not show low-energy quasiparti- 
cles (they exhibit exponential temperature dependence 
of the London penetration depth and no residual linear 
term in the specific heat), meaning that all bands are 
fully gaped. Other materials, with weak impurity scat- 
tering, like LaFePCi2l!, KFeaAsJ^l, BaFe2[Asi_2;P^]^, 
shows behavior characteristic of a clean superconductor 
with line nodes of the gap parameter, namely, linear tem- 
perature dependence of the London penetration depth 
and square root magnetic field dependence of the ther- 
mal conductivity. Moreover, for BaFe2[Asi_a;Pa;]2 the 
presence of the node lines in one of the bands has been di- 
rectly confirmed by ARPES'^S., Note that the existence of 
accidental nodes does not necessarily contradict the over- 
all picture of the s± state - even though the order param- 
eter may change sign within one band, what matters most 
is the sign of its average value over the Fermi surface of 
this band. The third group of materials is formed by the 
compounds with strong scattering by dopants, such as 
Ba[Fei_2;Cox]2As2, which have rather large residual spe- 
cific heat at low temperature^^H^ quadratic dependence 
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of the London penetration depttP^, and typically show 
residual thermal conductivity^. These properties are 
not consistent with fully gapped s-wave order parame- 
ter. While this behavior can be interpreted as indication 
for the accidental nodes of the gap in combination with 
strong impurity scattering, it is also compatible with s± 
state where interband impurity scattering generates large 
number of subgap states and leads to finite density of 
states at the Fermi leveP^. 

The discovery of another iron-based superconducting 
family^^ ~ the chalcogenide 122's - cast some doubts on 
the s± state as an universal state for all compounds, since 
in these materials the hole bands around the F point 
are absent. Possible pairing mechanisms and structures 
of the order parameter in view of recent experimental 
developments have been extensively discussed in recent 
review^. 




FIG. 1. On the left: The proximity structure composed of s- 
wave and s± superconductors which we consider in this paper. 
The superconductors are described by the Green's functions 
$ depending on coordinate x and Matsubara frequency u and 
corresponding gap parameters A (a;). For s± superconductor 
these parameters depend on the band index a = 1,2, see 
Sec. |III| for details. "l>s couples to the $q on the s± side. In 
the case of close coupling strengths with different bands, this 
creates the possibihty for a non-trivial relative phases (f)a (as 
functions of a; and uj), and the TRSB state; On the right: The 
three possible states. From top to bottom - the behavior of 
the anomalous part of the Greens' functions in the symmetric 
and asymmetric TRSB state, and the aligned state. 

In this paper we propose an alternative approach for 
probing the superconducting order parameter in iron- 
based materials. It is based on proximity effect - the mu- 
tual influence of two superconductors, brought in contact. 
We study the vicinity of a boundary between ordinary s- 
wave and s± superconductors, see Fig. [T] Competing 
interactions between the s-wave order parameter and s± 
gaps with opposite signs lead to frustration. As a conse- 
quence of this frustration, several interesting effects can 
appear, including the possibility of new superconducting 



states. This is in contrast with proximity between con- 
ventional (whether single- or multi-band) superconduc- 
tors, which is a rather straightforward phenomenon - the 
phases of the gaps on both sides always align and their 
amplitudes get closer (i. e., the smaller is enhanced and 
the larger is suppressed - see, for example, Ref l32|) . In 
the case of a contact between s and s± superconductors 
there is no obvious way to align the phases and several 
possibilities compete. 

Various effects close to a boundary between s and s± 
superconductors have been already considered. Both 
phenomenological and microscopic methods were used, 
number of interesting and novel results were obtained 
and different effects have been suggested as possible fin- 
gerprints of the s± state. Several papers have concen- 
trated on the problem of a Josephson junctions between 
s and s± superconductor^Ii^IMEll context of 

SQUIEP^ and JosephsorP ^*^ in terferometry, and macro- 
scopic quantum tunnelinj^lllll (for a detailed review of 
the theoretical and experimental results see Ref . f55|) . In 
Refs. [37H401 the mutual effects of s and s± superconduc- 
tors in contact were considered, and the possibility for a 
new time-reversal symmetry breaking (TRSB) state close 
to the interface was discussed. The methods used to ob- 
tain the TRSB state were different, however, all of them 
have some intrinsic limitations and the obtained solu- 
tions cannot be regarded as fully microscopic and self- 
consistent. 

In this paper we consider in detail anomalous features 
of system composed of a two-band s± superconductor 
in contact with a weaker (i.e., with smaller bulk criti- 
cal temperature) single-band s-wave superconductor, see 
Fig. [T] We employ simple microscopic model, which as- 
sumes dirty limit in all superconductors, but neglects in- 
terband scattering in s± superconductors. This model 
allows us to describe the system using a minimum num- 
ber of the most essential and physically transparent pa- 
rameters. Our approach neglects many features of iron 
pnictides, such as the presence of more than two gaps, 
interband scattering, orbital content of the bands, ac- 
cidental nodes etc., which may be relevant for detailed 
comparison of the theory with the experimental data. 
We believe, however, that the overall picture we present 
will survive even in a more realistic calculations. 

As we show below, several non-trivial effects can be ex- 
pected is such structures. If the s-wave superconductor 
is much stronger coupled to one of the gaps on the s± 
side than the other, it is natural to expect that its order 
parameter aligns with this gap. The phase of the gap 
on the other band is then anti-aligned with the s-wave 
gap and thus frustrated. This frustration leads to the 
possibility of negative proximity effect - the magnitudes 
of all gaps being suppressed close to the boundary (see 
the lower right panel on Fig. [l]). This effect IS unique 
to the interface with s± order parameter, and has to be 
contrasted with the case of s-s++-wave structures, where 
the smallest gap is always enhanced (we call this con- 
ventional proximity effect positive). As we demonstrate 
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below, in the aligned state both positive and negative 
effects are possible. 

When the properties of the two bands on the s± 
side are roughly the same, and the coupling across the 
boun dary is identical, even more interesting possibility 
arisea^ZESl Since the frustration in such system is large, 
another state emerges as a compromise. In it the relative 
phase of the gaps on the s± side deviates from its bulk tt 
value close to the boundary. This tilting relieves some 
of the frustration due to the inter-boundary coupling 
and disappears in the bulk (see the upper right panel 
on Fig. [T]) . Such state is intrinsically complex and thus 
breaks the time-reversal symmetry^. Due to the finite 
phase difference between the gaps, it is characterized by 
spontaneous supercurrents with opposite direction flow- 
ing between the s-superconductor and s± bands so that 
in ground state the total supercurrent is zero. In a more 
general situation of different s± gaps, the TRSB state is 
possible when the Josephson coupling energies between 
the s order parameter and the opposite-sign bands ex- 
actly compensate each other. It was also demonstrated 
recentljffl that a different kind of the TRSB state may 
realize near the free surface of s± superconductor with 
very strong interband reflection. 

Let us outline the structure of the paper. First, in 
section |ll] we present a very simple phenomenological 
model of a frustrated Josephson junction. In spite of 
its simplicity, we believe that this model catches some of 
the essential physics of the system. The results we ob- 
tain agree with the intuitive picture we presented above: 
the TRSB state is stable when the Josephson couplings 
with opposite-sign bands are very close corresponding to 
strong frustration. Away from this region, the aligned 
states are more favorable. The frustrated Josephson 
junction model predicts continuous phase transitions be- 
tween the aligned and TRSB states. 

The above simple model, however, does not describe 
proximity effects, and misses other important details as 
well. To develop a more realistic description of the super- 
conductivity on both sides of the interface, we use micro- 



from general considerations). In section VIII we summa- 



scopic theory in the dirty limit presented in sections III 



which describes the systems by multiband Usadel equa- 
tions supplemented with the appropriate boundary con- 
ditions. In section IV we present analytical results ob- 
tained in the limit of weak coupling between the super- 
conductors, and in section [V] the procedure we used for 
numerical solution of the equations is described. Even 
though this approach is strictly applicable only in the 
dirty limit, we expect our results to be qualitatively (or 
even quantitatively) correct in the clean case as well. 

Within the developed framework we study different 
superconducting states. We start with the more con- 



ventional aligned state (in section VI) and obtain both 
positive and negative proximity effects, depending on the 
values of different physical parameters. In section VII we 



show that the TRSB state indeed exists, and develop its 
quantitative description. Furthermore, the proximity ef- 
fect for such state is always negative (this is also clear 



rize our results and discuss the possible limitations of our 
approach. 

We again emphasize that the negative proximity effect 
is unique consequence of the bulk s± superconductivity 
- it cannot be present for conventional s-wave order pa- 
rameter (assuming that the interface does not have ex- 
trinsic pair-breaking surface effects). Thus, observation 
of such effect will be a definitive prove of the presence of 
s± state in iron pnictides. The reverse is not true - posi- 
tive proximity effect is possible for both s± and s_|__|_-wave 
superconductors. 

Another interesting consequence of our results is the 
possible phase transitions between these different super- 
conducting states. To induce the TRSB state we can 
vary the asymmetry. More relevant experimentally, how- 
ever, are the transitions tuned by the temperature. We 
expect the aligned state to be stable at higher tempera- 
ture, but be supplanted by the TRSB state (if it is the 
ground state for T = 0) at lower temperatures through a 
phase transition. 



II. SIMPLE MODEL: FRUSTRATED 
JOSEPHSON JUNCTION 

To gain some insight into the phase diagram of the 
system let us first consider a simple phenomenological 
model, which nevertheless catches essential physics - 
Josephson junction between s and s± superconductors 
(see also Ref. 130]). The energy of such systems depends 
on the phase shifts between the order parameters in the 
bands of the s± superconductor and the order parameter 
in the s-wave superconductor, 9i 2 = '/>i,2 ~ On gen- 
eral grounds, the simplest form of this energy in reduced 
form can be written as 

E{9i , 6*2) = cos(6'i - 6*2) - ti cos 0i - t2 cos 02 (1) 



with 



tn 



where Ej^a are Josephson coupling energies between the 
s-wave superconductor and two bands of the s± super- 
conductors, f 12 is the interband coupling energy, and d± 
is the thickness of s± superconductor. For definiteness, 
we assume ti < t2- Minimizing the energy, we get the 
equilibrium conditions: 

-sin(6'i - 02) +iisin6'i = 0, 
sin(6'i - 02) +^2 sin 6*2 = 0, 

and, as a consequence, 

ti sin 6i — ~t2 sin 62 ■ 

Excluding 62, we obtain equation for Oi. 



ti 



1 — ^ sin^ 61 — cos 01 — + 

t2 



= 0. 
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FIG, 2, Left part: The phase diagram of a frustrated Joseph- 
son junction. The TRSB state exists near the diagonal part 
where t\ ^ t2 ■ Away from this diagonal region the aligned 
state is realized. The right plot illustrates evolution of the 
phase angles 61^2 along the vertical line marked in the phase 
diagram. 



This equation has two solutions, corresponding to the 
aligned and TRSB states. For the aligned state, sin 9i = 
and we obtain 61 = Q, 62 = and 



4 -ii +t2 



For the TRSB state we obtain 
cos 01 — 

cos 02 = 



t2 
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t2 


2 " 


^ 2t^ 


2tl 


h 


1 


h 


2 " 


^ 2t[ 


24 



(2) 

(3a) 
(3b) 



In particular, for the symmetric case ti = ^2, and we have 
cos 01 — cos 6*2 — ti/2. The energy for the frustrated 
state can be written as: 



Eh = -^{tih+t2/h+h/t2). 



(4) 



The transition to the aligned state occurs when cos 61 — 
1, giving 



t2 



1 + h 



(5) 



and the ahgned state is stable for t2 < ii/(l + ii). Con- 
verting to real units, we conclude that for weak Joseph- 
son coupling the frustrated state is realized in the region 
\Ej,,-Ej^2\<Ej/{£ud±). ^ 

These results are summarized on Fig, [2] As we have 
anticipated, the TRSB state exists in the region where 
the frustration is maximum (ti w ^2), and is replaced 
by the aligned state when one of the couplings domi- 
nates. The transition between the two states is continu- 
ous. Even though this simple model can not pretend to 
quantitative description of the interface, we expect that 
its general features survive in a more microscopic setup, 
which we now proceed to describe. 



III. USADEL EQUATIONS AND BOUNDARY 
CONDITIONS 



We now move to a fully microscopic description of the 
problem. Let us consider a "sandwich" , consisting of a 
slab of two-band s± superconductor with thickness d±, 
in contact with slab of a single-band s-wave supercon- 
ductor with thickness ds , as shown on Fig, [T] We choose 
the a; = plane to be the boundary between the two 
layers. The main assumption of our description is that 
both superconductors are in the dirty limit, but the in- 
terband scattering in the s± superconductor is negligible. 
In this case superconductivity can be described by a sim- 
plified version of the Gor'kov equations, known as Usadel 
equation^^. The following formalism has been developed 
in Ref. and applied to the case of conventional two- 
band superconductors in Ref. 

For the s-wave superconductor {—ds < x < 0), the 
equations for the impurity averaged Greens' functions Gs 
and (where <I>s = ujEg/Gs and F, is the anomalous 
part of the single-particle Green's function) , with the nec- 
essary self-consistency equation, are: 



2wG< 



27rr VA,^^ = A„ G« 

^ ^ M 



LJ>0 



(6a) 
(6b) 



where uj = 2TiT{n + 1/2) stands for the Matsubara fre- 
quencies. 

For the s±-superconductor, Q < x < d± and the band 
index a = 1, 2, we have 



2ujGa 

2nT ^«/3 



A r' 



13, tJ>0 



V^2 + |$„|2 



(7a) 
(7b) 



(we again emphasize the fact that interband impurity 
scattering has been neglected - for details see Ref. ,3^ 
and Section Villi. We denote the bulk critical temper- 



atures of the s± and s-wave superconductors as and 
Tg*, respectively. Since we consider s± superconductor, 
we assume A1A2 < which is realized if Ai2,A2i < 0. 
The diffusion coefficients £'{s,a} are related to the con- 
ductivity a{s^a} as CJ{s^a} = e^V{8,a}D{s,a}, whcre V{s,a} 

is the normal density of states (DoS). The ratio of the off- 
diagonal coupling constants is given by the ratio of partial 
normal DoS, Xapl^pa = v^jva- It is convenient to nor- 
malize all energy scales (w and gaps on both sides) to ttTc. 
We also introduce coherence lengths — \J Da/2'iTTc 
and ^* = yjDs/2'n:Tc (note that ^* is related to the true 
bulk coherence length of the s-wave superconductor by 

Since we consider an interface, these equations have to 
be supplemented with appropriate boundary conditions. 
These connect the Green's functions and their derivatives 
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at the X = plane and can be written a; 



7Ba 



for a — 1,2 with 

7a 



7 7-Ba 7 ! 



(8a) 
(8b) 

(9) 



where Rsa is the boundary resistivity for band a. We 
can combine Eqs. (8a)-(8b) and get the useful equivalent 



form of the boundary condition 



(10) 



where we introduced the new interface parameters jBa — 
iBala which we will use together with 7bq. We also have 
to specify the conditions on the external boundaries: 



<i>;(-rf,) = 0, <(d±)-0. 



(11) 



Four parameters enter the boundary conditions and 
control the strength and the sign of the proximity effect - 
7i and 72 depend on the bulk properties of the materials, 
whereas 7^1 and 732 describe the boundary itselP^. The 
first two parameters determine the relative strength of 
the proximity effect between the s-wave superconductor 
and the s± bands. In particular, large 7^ implies that 
the s-wave material is more metallic than the a band 
on the s± side, and thus strongly influences it through 
proximity, while remaining weakly affected by this band 
itself. For the ratio of these parameters we derive the 
following relation 



71 
72 



^2^2 Ai2^2 



(12) 



The parameter 75^ is inversely proportional to the trans- 
parency of the boundary for the a band. Estimating 
these parameters is not easy, but for the case of iron- 
based materials (which typically are semi-metals) in con- 
tact with conventional superconductor we generally ex- 
pect 7q, to be large. 

To find the density of states, we have to perform ana- 
lytical continuation of the Green's functions to real ener- 
gies iuj ^ E + iS. The normalized DoS is related to the 
real-energy Green's function by the standard expression 



Re 



E 



E^-^{,^^y{E,x)^l^^{-E,x) 



(13) 



In the following section we present analytical results for 
the Green's functions and gap parameters obtained in the 
limit of weak coupling between superconductors (large 
7_Ba) in the case of aligned state. 



IV. ANALYTICAL RESULTS FOR WEAK 
COUPLING 

In the case of the weak coupling between s and s± 
superconductors, 7^0 ^ 1, contact-induced corrections 
to the gaps and Green's function can be treated as 



small perturbations, A{s q,} (x) = A 



{s,q}0 



A{s,a}(a;) 



${5.0,} (x) — A{s„}o + ^{s,a}ix)- Small corrections 
^{s,q}(w,x) and A^s a}{x) can be computed analytically 
in the linear order with respect to l/"fBa- Similar cal- 
culation for several types of junctions using somewhat 
different approach has been done in Ref. [371 We con- 
sider here the case of aligned gaps. The computation 
details are presented in Appendix [A| and general results 
can be presented in the form of Fourier expansions. For s- 
wave superconductor $s(a;,a;) = X]m=o^s,m(w) COS kfjiX 
and As(x) — X]m=o ^s,m cos kmX with fc„j — mir/ds and 
Fourier components 



■ 1+C 

E 



2 U2 



Afo A,,o 



f 2 1.2 

AqO 



2nT 



(14a) 



E 



(2 



5, 



^l^/ {d,Cs) Aso~ Am 



1 



iBa 



(14b) 



2 



-A^ 



3/2 



1 I t2 jU2 



where = Ds/{2^u^ + A%). 

For the s± superconductor the corresponding ex- 
pansions are = J2m=o^a,mCOsqmX, A^ = 
J2m=o^a,'mCOsqmX with Qm = rmT/d±. The Fourier 
components '^a,m and A^ „i are given by somewhat cum- 
bersome but closed analytical formulas, 

, (2-^™)e./(rf±ea) 



A. 



1 + £2 ff^ 



A^o A,o 



1 + 
Aqo 



2 „2 



+ A^o iBc 



,al3- 



(15a) 



,9,w>0 



A2 



(2-UClJ(rf±£/3)A,o 



1 



jBfS 



(15b) 



and £,1,^ = Dal{2^uj'^ + A^g). Here the matrix C/,„,„/3 
is defined by the following relations 
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W22-S™,2 -1012 

-W21 Wii-Sm^iJ 

2Wll-Sm,lti'22 + Em.l E77), 2 
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with 



Frustrated state 



LJ>0 



2 + A2„)^/' W 



-In 



and w — X ^ — X being the degenerate matrix whose 
components are given by relations 



Wll = 



W22 = 



A^/4 + Ai2A2i - A_/2 
dctA ' 
^A^/4 + Ai2A2i+A_/2 
dctA ' 



W12 = 



A12 
detA' 
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FIG. 3. The regions of positive and negative proximity in the 
weak-coupling regime evaluated using Eq. (16 1. The left and 



A12A2 



where A_ = An — A22 and det A = A11A22 

The quantity represents the average correction to 
the s-wave gap parameter induced by the contact. In 
particular, the sign of A^^o determines wether s-wave su- 
perconductivity is enhanced or suppressed (positive vs 
negative proximity effect) . At low temperatures it is pos- 
sible to obtain an analytical result for the average gap 
correction, see Appendix [A] 



right diagrams illustrate the cases of identical and different 
s± gaps. 

To find corrections to the densities of states, SNs{E) 
and SNa{E), we have to perform the analytical continu- 
ation iuj ^ E + iS in Eqs. (14a) and (15a) and perform 



expansion in Eq. ( 13 ) which gives 



with U{a) 



C 



u 



|Aao|/ 7Bq|Aqo| 

-E{\-a^) 



5N{,^a)o{E,x) = Re 



(16) 



{E-^ 



Ms, a} 



3/2 



(18) 



1 



where K{m) = J^^^{1 - msin^ e)-^/^d9 and E{m) 



Simple analytical result illustrating general trends can 
be obtained at low temperatures in the case of thin s- 
layer, dg < ^s, and weak s-wave superconductor, |Aq| ^ 
|Asp2l. In this case, the real-energy Green's function can 
be approximately evaluated as 



/J^^^(l — TO sin^ O^/'^dO are the complete elliptic integrals. $s(i? x) w A^ q — t:T^ 



1 



A 



This general results further simplifies for the important 
particular case A^q ^ |Aqo|- In this limit the elliptic 
integrals can be expressed in terms of elemental functions 
leading to 



(19) 



ttTc ds ^ 7sa|Aao| 



In 



4|A 



qOI 



A 



sO 



- 1 



Neglecting the first trivial term, we obtain correction to 
the s-wave DoS 



(17) 



2 n3/2 



From this result we can see that the partial contribution 
from the band is mostly determined by the strength of 
coupling to this band oc In addition, we observe 

that while the positive contribution from the aligned 
band is proportional to the gap difference Aio — A^g, 
the negative contribution from the anti-aligned band 
is proportional to the sum of the absolute gap values 
IA20I + AsQ. As a consequence, even in the aligned state 
the negative contribution may exceed the positive one 
leading to the total negative proximity effect. This the 
negative-proximity region is especially broad in the case 
when the gap values Aq,o and A^q are close. Figure [3] 
illustrates regions of positive and negative proximity in 
the coupling-constants plane, 1/7_bi - I/7B2, obtained 
using Eq. (16). The cases of identical and different s± 



{E^ - A 
^ Aqo — Aso 



e(|A,o|-£^). (20) 



where Q{x) is the step function. We see that the aligned 
bands (positive A^) induce positive corrections to the 
s-wave DoS and the anti-aligned bands (negative A^) 
induces negative corrections. These ne gat ive features can 
serve as definite fingerprint of s± state^. Note that the 



gaps are illustrated. We can see that in both cases the 
region of negative proximity occupies significant region 
in the parameter space. 



perturbative result ( 20 ) does not describes energy regions 
in the vicinity of the gap values, E ^ A^q- 

To go beyond the weak-coupling regime we have to 
rely on numerical calculations. In the following section 
we will describe numerical procedure and present results 
of these calculations for different cases. 



V. NUMERICAL PROCEDURE 

For numerical modeling, it is convenient to use the so- 
called 0-parametrization, in which we write Gg = cos 9s , 
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Ga = COS Oa- For $Q and $s two choices will prove conve- 
nient. Let us first consider the (technically simpler) case 
of a significant difference between the coupling of one of 
the two gaps on the s± side to the gap on the s-wave 
side (weakly frustrated interface). As already mentioned 
in the introduction, in such situation we expect the so- 
called aligned state to be stable. In it the phase of $s 
is aligned with one of while the phase difference be- 
tween $1 and $2 is TT - this state obviously belongs to the 
general class of s± states. Proximity effect on the s-wave 
side enters through the (asymmetric) suppression of Ai 
and A2 close to the boundary (schematically shown on 
the lower right panel of Fig. [T]) . 

In this state it is possible to choose $q, and $5 real, 
^{s.cx} = w tan 6'{3 ct} (with the condition <I'i$2 < 0). In 
the following, we assume that the coherence lengths of 
the s± bands are equal ^1 — £_2- We use dimension- 
less units: all energies are normalized to ttTc and the 
lengths on right/left are normalized to /^q. Separating 
the highest derivatives 9'! , , we can rewrite the Usadel 
equations as 



^{s,a} 



+ A 



{s,a}' 



7{s,q} 



io sin0 



0, 



(21) 



This equation determines 0{s.a} as function of discrete 
Matsubara frequency w = t{2n + 1) with t = T/T^- 

For simplicity, we neglect the intraband pairing inter- 
actions and consider only the repulsive interband cou- 
pling, parametrized by A12 and A21. This is a reasonable 
simplifications for the case of iron pnictides, in which the 
interband pair-scattering is believed to be driving the 
superconductivity. Proper generalization, including the 
(attractive or repulsive) intraband terms, is straightfor- 
ward. The self-consistency equation for A^ ( 7b ) becomes 



Aq, = —2t\af} smOjj. 



(22) 



n=0 



with a, /3 = 1, 2 or 2, 1, iVn,ax = a;max/(27rr) - 1/2, where 
Wmax is some frequency cutoff. 

The boundary conditions at the outside boundaries, 
X = —ds,d±, are 9'g{—ds) = and 9'^{d±) = 0. For the 
boundary conditions for the aligned state at the s-s± 
interface we obtain (see Appendix IbI) 



E 



la 



sm 



iBc 



(23a) 
(23b) 



Let us now consider the case of strongly frustrated 
boundary when the TRSB state appears. The Green's 
functions are now essentially complex quantities - '^s — 
uj tan 9 sc'^'^" and = tan 6*06"''° (the factor i in the 
definition of is for convenience). For the s± state 
we have ipi — — tt. Correspondingly, the s-wave gap 
and the two gaps of the order parameter can be writ- 
ten as AgC^^' and iAae^^"" , where the (real) parameters 



and X{s,a} have to be determined from the self- 
consistency equations. 

The transformed Usadel equations for the highest 
9" 



derivatives 9'! 1 and lo'! , in this case become: 



(p" + 2 cot ( 



- sin((/3 - ^) = 0, 



(24a) 



sin 9 

9" — sin 9 cos 9{(p'f-ujsm9 + A cos 9 cos((/7 — x) = Oi 

(24b) 

where for brevity we omitted the subscripts {s,ck}. We 
use the same normalization for all energies and lengths 
as in Eq. pTl). 



The self-consistency equation for Aq is: 



AqC*-^" = —2t\ap ^ sin^^e 



(25) 



This complex equation can be split into two real equa- 
tions for Aq and Xa'- 



Aa = 

2t\al3 

Xa = - 



sin 9p cos (^/3 I + ( sin 9p sin ipp ] 



- arctan j sin 9p sin tpp / sin 9p cos Lpp 



1/2 



The negative sign in the second equation - a straight- 
forward consequence of the fact that Aq^ is repulsive - 
leads to somewhat counterintuitive result. The proxim- 
ity effect tends to align the phases of $s and 
in turn, determines the order parameter A^ through the 
self-consistency equation, which tends to make the phases 
difference between $q, and A^ close to tt (to compensate 
for the "wrong" sign of the interaction) . This means that, 
due to the proximity effect, the phase difference between 
Aq, and As increases (measured from the positive axis), 
rather that decreases - see Fig. [l] for a pictorial presenta- 
tion of this argument. This is a somewhat oversimplified 
picture, since </? depends on Matsubara frequency. As 
can be seen from the Eqs. (7a)-(7b), for high frequencies 
the derivative term becomes unimportant and Aq, 
and (pa — ^ Xa- This implies that (fa changes its sign 
as a function of a;. To obtain A/3 we have to sum over 
all n and this leads to "smearing" of the contributions 
to xi3 over a phase interval. Nonetheless, as we will see, 
the intuitive argument above seems to be qualitatively 
correct. 

The boundary conditions at the outside boundaries are 
given by 

0',(-4) = O, p's{-ds) = Q 

The boundary conditions at the s-s± interface in 9- 
parametrization are more complicated and we derive 
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them in Appendix [B] In reduced units these conditions 
are 

1 



iBc 



[cos 9aSm0sSn\{ips — 'fa) — COS 9s sin 9 a] 



Lp'^ sin 9a — sin 9^ cos ((^^ — ipa) ■ 

iBa 



(26a) 
(26b) 



9's — — — [cos 9a sin 6's+cos 9s sin^jj sin [ipa — ^s)] , 

(26c) 
(26d) 



(f's sin 9s — sin 9^ cos {ipa — (fis 



To summarize, we have developed an entirely self- 
consistent scheme, based on Usadel equations and sup- 
plemented by the appropriate self-consistency equations 
and boundary conditions. Unfortunately, due to its con- 
siderable complexity, it has to be solved numerically in 
the general case. To do this we start with a guess for 
9s, 9 and if, solve it on the s± side, using two of the 
boundary conditions at a; = and the boundary condi- 
tions at X = d±, then write the third condition at x = 
with the new solutions for 9, ip and use it to solve for 9s. 
With such obtained solution we re- write the two initially 
used boundary conditions and again solve for 9a and (pa ■ 
We repeat the process until self-consistency is achieved. 
To simplify the calculations, in some cases we expanded 
the equations around the bulk value 9a = arctan(Ao/w), 
ip = 0. 

Once calculations have produced self-consistent solu- 
tions of the equations on both sides of the boundary, we 
can obtain the DoS at T = 0, N^s,a}{x)- For this, we 
rewrite the equations and boundary conditions with an- 
alytically continued frequency iiv ^ E + id, and with 
computed gap functions As{x), Aq,(x) and solve them 
again for the real-energy Green's functions 9^s.a}{E,x). 
Once these equations are solved, the DoSs given by Eq. 
( 13 1 or, in the 0-parametrization, by 



N{s,a}iE,x)=Re[cos 9{s,a}iE,x)]. 

In the following sections we well review properties of the 
proximity system for different cases. 



VI. ALIGNED STATE 

With the formalism described in the sections |III[ |IV| 
and |V] we are now ready to study the proximity effects 
in several particular cases. 

A. Identical s± bands 

We consider an asymmetric s± state. By this we 
mean that the properties of the bands and the bound- 
ary condition parameters are not symmetric with respect 
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FIG. 4. These plots illustrate behavior of the gaps in the 
proximity sandwich for two values of the ratio 7S1/7B2. The 
insets show blowups of the As [x) and for comparison A^""^ 
is also shown. Other used parameters are 7si ~ 10, ds = Csi 
and 7a — 10. One can see that for weakly asymetric coupling 
7si/7s2 = 0.75 in the plot (a) the proximity effect is neg- 
ative (As (a;) < As""^) while for strongly asymetric coupling 
7si/7s2 = 0.25 in the plot (b) the proximity effect becomes 
positive (As(x) > A^""'). 



to the band indices interchange. For simplicity, we first 
consider the asymmetry only in the coupling strength 
7-Bi 7^ lB2, while keeping the rest of the parameters 
symmetric. In particular, this means that close to the 
boundary |Ai| ^ IA2I, but in the bulk the symmetric s± 
state is restored and A2 — — Ai (see right middle panel 
of Fig. [I]) . In all numerical calculations presented in this 
section, unless stated otherwise, we fixed several param- 
eters: d± = 8^a, = 0.222Tc, T = O.lTc. The choice 
of gives relation for the gap A^""^/(7rrc) = 0.1213. 
Other parameters which enter the calculations are spec- 
ified in the figure captions. 

Since in a way the aligned state is closer to a conven- 
tional multiband superconductor than the TRSB state 
is, we expect the usual positive proximity effect to be 
predominant. Indeed, as we will see, the s-wave super- 
conductor's gap tends to be enhanced by the presence 
of the stronger s± superconductor (which in turn is sup- 
pressed). This is easy to understand when the s-wave 
gap is much stronger coupled to one of the gaps on the 
s± side, since this case can be thought as proximity be- 
tween two single-band superconductors. Coupling to the 
other gap on the s± side can be treated as a small per- 
turbation. However, when the interboundary couplings 
are close, and the system is strongly frustrated (and thus 
close to the TRSB state) , the proximity effect turns neg- 
ative, with all superconducting gaps suppressed close to 
the a; = plane. Thus the "sign" of the proximity effect 
on the s -wave side in the aligned state is not universal, 
and depends on both bulk and boundary properties of 
the compounds. 

On Fig. |4] we illustrate this effect by showing calcu- 
lations of As{x) for two different values of 751/732 (we 
change 7^2 while keeping all other parameters of the sys- 
tem fixed). As can be seen, A^ is enhanced or suppressed 
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FIG. 5. (a)Plot of As(— ds) and As(0)) as a function of the 
ratio ^Bil^Bi- As ^b\/^b2 approaches 1, the enhancement 
on the s-wave side is replaced by suppression (the proximity 
effect changes from positive to negative). The arrow indicate 
the expected location of the transition to the negative proxim- 
ity from the weak-coupling approximation, 731/7^2 ~ 0.645. 
(b,c) The dependences of the s-wave gap parameter at the 
boundaries (plot (b)) and the a± gaps at the interface (plot 

(c) ) on the coupling strength with the first band, I/7S1 for 
fixed ratio 7_b2/7si = 2. The inset in the plot (b) blows 
up the small-coupling region. Behavior of the gaps expected 
from the linear approximation is also shown (the dotted line 
in the plot (b) and the dashed and dotted lines in the plot (c)). 

(d) The regions of the negative and positive proximity in the 
I/7S1-I/7S2 plane obtained by numerical calculations. The 
dashed lines show boundaries obtained within weak-coupling 
approximation. Parameters used in the calculation for all 
plots are 7^1 = 10, 7^ = 10, ds — Q- 



close the interface (corresponding to a sign change of 
6*^(0)), depending on the ratio of the boundary trans- 
parencies. Note that the transition between the two cases 
coincides with change of the ratio of Ai(0)/|A2(0)|. For 
positive (negative) proximity effect this ratio is larger 
(smaller) than 1. This is easy to understand from the first 
condition in Eqs. (23a) - it is clear that (for 71 — 72) the 
sign of O'g is determined by the 6'i + 9'2- Also remember 



that Ai is determined by 62 and vice versa. 

Now we proceed to systematically study the interplay 
between the different physical parameters and the tran- 
sition from negative to positive proximity effects. We 
plot on Fig. [sjja) the value of As{—ds) and As(0) as a 
fimction of the ratio 7bi/7B2 (the model, of course, is 
symmetric with respect to the exchange 7^1 O 752 )• 
As is suppressed as this ratio gets closer to one (and 
the frustration increases), in agreement with our previ- 
ous qualitative arguments and explicit calculations. With 
the decrease of the ratio the proximity effect turns from 
positive to negative when the ratio 751/732 exceeds the 
critical value sa 0.578. The linear approximation for the 



average correction to the s-wave gap ( 16 1 in the case 
IA20I = Aio gives the following estimate for this criti- 
cal value, 7Bi/7s2 ~ (Aio - Aso)/(Aio + A^o). For our 
parameters this gives 7bi/7B2 ~ 0.645 which somewhat 
exceeds the value obtained in numerical calculations. As 
the system gets closer to the line of maximal frustration 
7bi = 7b2 it eventually undergoes a phase transition to 
the TRSB state. 

On Figs. [5jb) and [5|c) we illustrate dependence of 
the gap parameters on the coupling strength 1/751 for 
fixed ratio 752/751 = 2. At small coupling strength, 
1/751 < 0.07 behavior of all gaps agrees with the linear 
approximation described in Sec. IV At larger coupling 



the s-wave gap strongly deviates from the linear approxi- 
mation. In particular, it reaches maximum at some value 
of coupling and decreases at larger values. With fur- 
ther increase of coupling, Ag drops below the bulk value, 
i.e. proximity becomes negative. Qualitatively, this be- 
havior can be understood as follows: The positive and 
negative contributions to A^ are roughly proportional to 
Ai(0) — As and |A2(0)| -|- A^ correspondingly. As s± 
gap parameters reduce and become closer to the s-wave 
gap with increasing coupling, the negative term becomes 
relatively stronger. Another notable property is that the 
absolute values of the s± gap parameters Ai and IA2I 
remain very close, in spite of significant asymmetry in 
coupling strength. 

On Fig. [5]jd) we show the regions of the negative 
and positive proximity in the coupling-strengths plane 
obtained by numerical calculations and compare them 
with predictions of the weak-coupling approximation. 
In agreement with the plots [sja-c), we can see that 
stronger coupling favors negative proximity meaning that 
the weak-coupling approach underestimates the width of 
the negative-proximity region for small 7bq. 

Unfortunately, there is no obvious way to control the 
coupling parameters 75Q, in experiment. One parameter 
which can be varied relatively easy in experiments is the 
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FIG. 7. (a)The coordinate dependences of the s-wave order 
parameter for different couplings to tlie s±-bands. The 
plots are labeled by 751,732- Other parameters are listed 
in the plot (c). The horizontal dashed line shows bulk gap. 
The curves above this line (positive proximity) correspond to 
the alignment with the large gap, Ai, while the curves below 
this line (negative proximity) correspond to the alignment 
with the small gap, A2. Three dot-dashed line above and 
three dashed lines below show ana lytical predictions based 
on linear approximation, Eq. ( |14b[ ). (b) The dependences 
of the s-wave gap at the outside boundary on the coupling 
strength with the aligned gap dashed lines show linear 

approximations, (c) The coordinate dependences of the s±- 
wave order parameters for the alignment with the larger gap. 
The solid lines show the linear approximation. 



thickness of the layers. On Fig. [6] we show the As as 
a function of dg for 7^1 = 10 and two values of 7b2j 12 
and 20, corresponding to negative and positive proximity 
effect. There are two natural tendencies which can be ob- 
served. First, with increase of the value of A^ increases 
(decreases) for the positive (negative) proximity. Second, 
the difference between As(0) and As{—ds) increases with 
thickness until it finally saturates for dg ^ ^s- 

In this section we demonstrated that both negative 
and positive proximity effects may present close to s-s± 
interfaces for the aligned state. The sign of the prox- 
imity effect is determined by non-trivial interplay of the 
physical parameters in the system. In general, negative 
proximity effect is expected in the region around the line 
which separates the aligned and the TRSB states, for 
intermediate frustration of the interface. 



B. Nonequal s± gaps 

In this section we present numerical results illustrating 
properties of aligned state when the bulk gaps of s± su- 
perconductor have different magnitudes. This situation 
is probably more typical and has one definite practical 
advantage with respect to the case of identical or close 
s± gaps - different gaps induce features into the s-wave 



DoS which are well separated in energy and thus easier to 
detect experimentally. Therefore, in this section, in ad- 
dition to behavior of the gaps, we also study behavior of 
DoS. We investigate in detail the shapes of DoS features 
and their sensitivity to the coupling parameters. 

We first consider the case of weak coupling (large 75^ 
and 7bq) and compare numerical calculations with the 
analytical results presented of Sec. IV In particular, this 
allows us to evaluate limits of the weak-coupling approx- 
imation. Fig. [7] illustrates behavior of the gaps for differ- 
ent coupling strengths. For s± superconductor we con- 
sider again the interband coupling model and the gap 
values are fixed by the coupling constants which we take 
as A12 = —0.4 and A21 = —0.2 giving the bulk gaps 
|Aio| = 0.683 and IA20I = 0.456 (in units of ttT^). We 
assume tcs — T^/Tc = 0.3 giving A^o — 0.1675. We again 
assume identical coherence lengths for two s±-bands. 
Due to relation ([l2]) we have 71/72 = A12/A21 = 2 and 
therefore we select 71 = 4, 72 = 2. All parameters used in 
these calculations are listed in the plot [tJJc) . We consider 
both cases of alignment of the s-wave gap A^ with large 
gap Ai (As tt Ai) and with small gap A2 (As TT A2). 
We only vary coupling strengths oc preserving the 

ratio 7'B2/7si = 2 for As TT Ai and 751/752 = 2 for 
AsTTA2. 

Figure [7|[a) shows the coordinate dependences of the 
s-wave order parameter. The symbols show numerical 
results and the lines show analytical results in the linear 
order with respect to l/7Ba based on Eq. (14b). The 
bulk gap is shown by the dashed line. For selected pa- 
rameters, the alignment with the large gap corresponds 
to positive proximity, while the alignment with the small 
gap corresponds to negative proximity. We can see that 
the linear approximation accurately describes behavior 
of As only for weakest coupling 7b(i,2) — fOO. It always 
overestimates the gap parameters and breaks down al- 
ready for rather weak coupling strength. This is even 
clearer on Fig. [Tjjb), where the gap parameter at the out- 
side surface is plotted versus the coupling strength with 
the aligned gap. Deviations are especially large in the 
case of alignment with the large gap. The found non- 
monotonic dependence of As is similar to the case of 
identical s± gaps shown in Fig. [5]Jb). 

Figure [tJ^c) shows the coordinate dependences of the 
s± gaps in the case As TT Ai. As expected, both gaps are 
suppressed at the interface. This suppression grows with 
increasing coupling strength, and disappears away from 
the boundary for a distance of the order of the coherence 
length. Behavior for the case As TT A2 is very similar. It 
is peculiar that, in contrast to the s-wave gap, the linear 
approximation accurately describes behavior of the s± 
gap parameters in the whole studied range of parameters. 

Figure|8]illustrates behavior of the DoS for the same set 
of parameters. The plots jsj^a) andjsjjb) demonstrate full 
shapes of the s-wave DoS, Ns{E), for the cases AsTTAi 
and As TT A2 correspondingly and for strongest coupling 
in this series. Interaction with the s± superconductor 
induces specific features in the s-wave DoS near the en- 
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FIG. 8. Behavior of the densities of states for the same pa- 
rameters as on Fig. [7] The plots (a) and (b) iUustrate shapes 
of s-wave DoS for the cases when the s-wave order parameter 
is ahgned with large and small gap correspondingly for the 
strongest used coupling strength. The dashed lines show the 
BCS density of states. The plots (c) and (d) show the par- 
tial DoS for two s± bands at the interface and at the outside 
boundary for these two cases. The plots (e) and (f) demon- 
strate the evolution of the feature in the s-wave DoS with 
varying strength of coupling. The plots are displaced verti- 
cally for clarity. The dashed lines show predictions of the 
linear approximation. 



ergies of s± gaps. Figures |8je,f) show evolution of these 
s-wave DoS features with increasing coupling strength. 
The weak-couphng approximation for the case of thin s- 
superconductor layer, Eq. (20), suggests that the aligned 



and anti-aligned gaps should induce a small peak and 
dip correspondingly^^. For comparison, the analytical 
weak-coupling results are also shown in the plots [8]je,f) 
by dashed lines. We can see that, the analytical approx- 
imation describes well the shapes of correction except 
the regions close to the s± gaps where it overestimates 
the amplitudes of peaks and dips. Nevertheless, small 
asymmetric peaks and dips do appear for weak coupling 
strength, in agreement with analytical predictions. These 
peaks and dips are rapidly smeared with increasing cou- 
pling strength. Eq. (20) also suggests that the correc- 



tion is strongly asymmetric: the s-wave DoS is only en- 
hanced or reduced for E < |Aq,|. This steplike behav- 
ior of the correction is indeed seen in the plots |8|^a,b) 
and, in more detail, in plots [8][e,f) . In fact, with increas- 
ing coupling the peaks and dips evolve into smooth up 
and down steps. The amplitude of this steplike feature 
between IA2I and |Ai| monotonically increases with in- 




FIG. 9. Evolution of the s-wave the order parameter (left 
lower plot) and DoS (left plot) for the fixed total boundary 
resistance set by condition I/7S1 + I/7S2 ~ 0.15 and for vary- 
ing partial resistances for two bands. Coupling to the larger 
gap in the right plot progressively increases from the back to 
front plot. Plots are marked by 7si, 732- Other parameters 
are the same as in the previous plots. The upper left plot 
shows the order parameter at the outside boundary and the 
gap in spectrum (defined in the right plot) as a function of 
a fraction of the first-band conductance with respect to the 
total conductance through the boundary. 



creasing coupling and its shape is well described by the 
analytical approximation. 

The plots |8jc) and [sj^d) show partial DoS for two s± 
bands at the interface and at the outside boundary. Cou- 
pling with s-wave superconductor leads to considerable 
smearing of the DoS peaks at the interface. The most 
prominent feature is the appearance of the tails spread- 
ing down to the s-wave spectral gap. This is a well-known 
feature which is always induced in a superconductor by 
proximity with either metal or weaker superconductopSB. 
Note also that the peaks positions for DoS at the inter- 
face do not shift much with respect the bulk gap values. 
There is no obvious difference in the DoS shapes between 
the aligned and anti-aligned bands. 

We now vary the partial resistances for two s± bands, 
while keeping fixed the total boundary resistance. As 
the total boundary conductance — ^/Rb is equal 
to sum of the partial boundary conductances l/Rsa oc 
I/tboi the total boundary resistance can be set by fix- 
ing 1/jB = I/7B1 + 1/7S2- We select this parameter 
as I/7B = 0.15 corresponding to the moderately strong 
coupling strength. All other parameters are the same as 
in Fig. [7] Figure [9] illustrates evolution of the s-wave or- 
der parameter and DoS for varying partial resistances for 
two bands. The plots are marked by the pair 7bii7b2- 
We can see that for strongly asymmetric coupling prox- 
imity effect is positive, there is pronounced BCS peak 
near the bulk gap, and pronounced steplike features at 
the s± gaps. The step amplitude is typically larger for 
smaller gap. As coupling becomes more symmetric, the 
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proximity turns negative, the BCS peak smears, and the 
amphtudes of the proximity-induced features decrease. 
Also, the gap in the spectrum corresponding to vanishing 
of DoS, Gs, rapidly decreases as coupling becomes more 
symmetric. This is seen more clearly in the upper left 
plot where order parameter and spectral gap plotted as 
function of a fraction of the first-band conductance with 
respect to the total conductance, / = Rb/Rbi- Both pa- 
rameters have quite sharp cusplike dependence on / and 
the spectral gap practically vanishes in the minimum. 
Note that, in contrast to the BCS DoS, the spectral gap 
is significantly smaller than the order parameter. 

In the next section we consider behavior of the order 
parameter for the TRSB state. 



VII. TRSB STATE 

For illustration we consider here only the simplest 
TRSB state, in the case the two-gap superconductor 
has identical bands, coupled symmetrically to the s-wave 
layer (see the panel at the right upper corner of Fig. [T]). 
The symmetry simplifies the calculations considerably. 
This is also the case which favors the TRSB state most, 
since the frustration at the interface is the strongest pos- 
sible. 

Due to the symmetry of the problem we can choose 
$3(0) to be real. It is also clear that $^ — $2 and 
we can write $1 — icj tan , $2 = ~iuj tan Oe^^"^. 
However, to preserve the uniformity in notation we keep 
the (now redundant) a indexing - (pi ^ —ip2 = (p, 
6*1 = —62 = 0. The above implies, of course, that we can 
write the gap functions as Ai = Ae*'<^, A2 = — Ae~*-^. 
Lp and X parametrizes the deviation from the s± state 
and we expect that in the bulk if ^ Q and $i,$2 be- 
come purely imaginary. From the boundary conditions 
it is clear that $5 can be chosen real everywhere (both 
lm($s(0)) — lm($'j,(0)) = - no imaginary component 
develops). In the 0-parametrization the s-wave Green's 
function is again determined by Eq. (21 K The equations 
for the s± side are identical to Eqs. 



(24 



but there are 

now only two independent variables (instead of four). In 



this case the boundary conditions ( 26 1 can be simplified 

as 



1 

1 



sm Us cos 6 a sm ipa + cos Os sm 6 a 



— sm Ws CSC f „ COS ipct 

iBa 

2 

6'g = COS 9s sin 9 a (tan 9s cot 0^ + sin ipa) ■ 

lalBa 

We proceed to solve the equations numerically. On 
10 we show the results of calculation for A, x 
again remind the reader that A2 — — Ai, X2 — —Xi) 
As done for different values oi^Ba- Several things should 
be noted: First, x(a;) is positive and thus the s± gaps are 
pushed away from A^, in accordance with the argument 



Fig. 



Aj2=|Ajexp[±/(7r/2+x)] 




FIG. 10. Behavior of the order parameters in the symmetric 
TRSB state for three different boundary transparency param- 
eters and for 7a = 5. The rest of the parameters are given 
in the text. Lower plots show the coordinate dependences of 
the absolute values of the order parameters and upper plot 
show the phase x{^) of the s± order parameter. Xa goes to 
zero in the bulk, where the state is s± , with phase rotated by 
7r/2 with respect to the s-wave OP. Note that jAi(0)| and Aa 
decrease and the value of Xa (0) increases with the increase of 
the interface transparency. 



in Section III Second, several away from the bound- 
ary, X ~ and the order parameter returns to its s± 
bulk form (since the TRSB state exists because the frus- 
tration created by the boundary); Third, decreasing ^Ba 
(increasing the boundary transparency) tends to reduce 
both Aq(0) and As(0), and increase x(0)i as expected. 
It is important to note that decreasing ^Ba pushes As(x) 
down everywhere, while keeping its general shape intact 
(unlike the changes of Aq(x)). This is due to two fac- 
tors, which generally decrease the value of A^(0), and 
thus enhance the strength of the proximity effect away 
from the interface. The first one is the relative thinness 
of the s-wave layer, and the other is the large value of 

7a- 

Since controlling the interface parameters ^Ba is not 
easy experimentally, we also study evolution with in- 
creasing thickness of s-wave superconductor ds for the 
same parameters of the interface. We vary the thickness 
of the s-wave layer from ^* to 4^* and show the results 



of the calculation on Fig. 11 On the s± side the changes 
are modest - x shifts a bit, while A^ is virtually un- 
changed for different dg- Note that as ds decreases, the 
gap on the right side gets slightly closer to the bulk s± 
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0.09 



-4 -3 -2 -1 

0.0 1 ' I 

-4-3-2.-1 1 2 3 4, 5 6 7 8 

FIG. 11. A(a;), x(a;) and As(a;) for = 5,7i3c, = 0.75 and 
four different values of ds- The changes on the right side are 
generally quite small. In contrast, As (a;) is suppressed strongly 
for thin films. This is more clearly seen in the lower inset 
which shows blowup plots of As{x). The upper inset shows 
the da dependence of the As at the outer surface. 

state (x goes down). The changes on the other side are 
more pronounced. is always suppressed as a: — ^ 0, 
but, as in the previous case, is significantly below its 
bulk value everywhere. 

These results are summarized in the upper inset of Fig. 
[TT] where we plot the value of the gap on the external left 
boundary as a function of the s-wave layer thickness. We 
see the result of the negative proximity effect as a sup- 
pression of As(— rfs) for thinner films. With the increase 
of the film thickness As(— dj) goes up, however, it stays 
noticeably below its bulk value even for dg = 4^* . This is 
again due to the relatively large 7q, which enhances the 
proximity effects away from the interface. 

Thus we can conclude that in the TRSB state the nega- 
tive proximity effect is present and, depending on the pre- 
cise values of various physical parameters, can be quite 
pronounced. 



VIII. DISCUSSION AND CONCLUSIONS 

One important point should be addressed before we 
apply our conclusions to structures involving iron-based 
superconductors. These materials have quite short co- 
herence length and are likely in the clean limit. In con- 
ventional superconductors there is smooth interpolation 
between the clean and the dirty limit and results ob- 
tained by solving Usadel equations are expected to be 



qualitatively or even quantitatively correct even in the 
clean limilpS'. Furthermore, the boundary itself acts as a 
scatterer, and can push the region up to several ^ away 
in the dirty regime. This again justifies the use of Usadel 
equations. In s± superconductors the situation is much 
less clear - because of the unconventional nature of the 
order parameter interband impurity scattering is pair- 
breaking. If strong, it can completely destabilize the s± 
stat^^. However, if this scattering channel can be ne- 
glected, as we have done (so s± is stable even close to 
the interface), then our approach should provide a rea- 
sonably good description. Whether this can be justified 
in realistic experimental situations is unclear, at least at 
the moment. 

In conclusion, we have studied the proximity effects 
close to a boundary between s and s± superconductors. 
Based on frustrated Josephson junction model, we have 
suggested phase diagram for such system. Because of 
the frustration, present at such interface, several inter- 
esting phenomena are possible. In the case of maxi- 
mum frustration - when the coupling of the two gaps 
on the s± side with the s-wave superconductor is compa- 
rable - a state which breaks time-reversal symmetry ap- 
pears. Such superconducting state is also characterized 
by a negative proximity effect - because of the frustra- 
tion the gap amplitudes on both sides are suppressed. In 
the case of significant asymmetry in the inter-boundary 
coupling s± can be stabilized even close to the inter- 
face, but becomes asymmetric there. Interestingly, such 
aligned state may also lead to negative proximity effect. 
It is very important to note that this effect is unique to 
the s-s± structures, and is not present close to a con- 
ventional S-S++ boundary. Observation of such effect 
in structures with iron pnictides/chalcogenides will be a 
definitive proof that their order parameter belongs to the 
unconventional s± class. 
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Appendix A: Boundary-induced corrections in the 
the weak-coupling limit for aligned state 

In this appendix we consider weak-coupling case 
iBa S> 1 and derive corrections to the Green's functions 
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and gap parameter for aligned state in the linear order 
with respect to the coupling parameters 



s-wave superconductor 



while the second contribution ^s,A obeys the following 
equation and the boundary condition 



(A5a) 
(A5b) 

The solution ^s,b{x) of the linear equation (A4a) with 
the boundary condition <1>^ = at a; ~ —ds is given by 



The first-order corrections to the s-wave Green's func- 
tions $s and gap As obey the following equations 



A, 



(Ala) 
(Alb) 



with 4% = D,/(2y^j^TAf) = 

and A ~ -^s/ (2As). In the boundary condition wc 
can neglect difference between <i>'s and A's and approxi- 
mate A's by their bulk values. This gives 



(A2) 



at X = with 7bq, = jajBa- The self-consistency condi- 
tion also can be rewritten as 



$s,&(a;) = Cs.hcosh 



X + ds 



(A6) 



where the constant Cs,b can be found from the boundary 
condition at a; = 



sinh(4/^s,^) ^ 7_bqG, 



E: 



Go 



-(As-A„). 



We compute ^s.A and Ag using the Fourier ex- 
pansiouj ^s,A = S^^o ^''.A,™ cos fc^x, A^ = 
' }2m =n As, TO COS kjnX with = rmr/ds. For ^s,b{x), Eq. 
iA6h, the Fourier components 'ts.b,m can be computed 
explicitly 



{2-S„,)es,J{dsCs) ^ yZJ^+Af As-A„ 



1 I f 2 )L2 



^ ^ (A7) 

Eq. (A5a| immediately gives the following relation be- 



tween the Fourier components $s.A,m and As 

As.m 



S , A , 771 



1 + ^l.uj^m 



(A8) 



cj>0 L 



(c.2+A2)3/2 



As 



As In ^ 0. 



(A3) 



We split 3>s hito the two contributions, $s = ^'s.fc + *fs,Ai 
where $s,6 is induced by the boundary condition and 
^'s,A is induced by the gap correction. The first contri- 
bution $s,b can be found from the following equation and 
the boundary condition 



lu^^'Lb - ^sM = 0, 



Ga 

iBa 



(A4a) 

(As - A„), (A4b) 



Subs titutin g this result into the self-consistency condi- 
tion (Alb I, we express As,m via ^s,b,m 



As,TO = {27tTJ2 



^0 {c^ + A^gf 



- s.b.m 



^„(^2 + A2)3/^' 



(A9) 



Eqs. (A7), (A8|, and (A9| provide a complete solution 



of the problem. 

For convenient comparison with the numerical calcu- 
lations, we also present these results in the reduced form 
in which dg is measured in units of = Dg/{2'KTc) and 
energies (w, A's and $'s) in units of ttTc 
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A, 



+ A2) 



3/2 



2^ E 



{-Km/ds) 

Vc<^2 + A2 + (7rm/4)' 

(2 - <5™)/4 A„ - A, 



l + (^m/4) /v^2 + a: 



2 - S,, 



^ . /7;2~ 



A„ - A, 



(AlO) 
(All) 



with t = T/Tc and w ~ 2t{n + l/2). At low temperatures 
the summation with respect to the Matsubara frequencies 
can be replaced by the integration ZttT^^^q — >■ dui. 
In this limit we can obtain an analytical resuh for the 
average correction to the order parameter, A^^Qj 



with U{a) 



A..0 ^ 



(A12) 



(z2 + 1)3/Vfl2z2 + 1 

X(l-a2)-j;(l-a2) 
l-a2 ' 

r7r/2 



where X(to) = (1 - m sin^ 6*) "1/2^^61 ^nd E{m) = 

/J^^^(l — TO sin^ OY^'^dO are the complete elliptic integrals. 

We present also the simple analytical results for im- 
portant particular case of (i) thin s-layer, ds ^ S,s,i\, 
(ii) weaker s -superconductor, Aj ^ |Aq|, and (iii) low 
temperatures, T ^ T^. Due to the first condition, the 
dominating contribution to the gap correction is given 
by the coordinate independent part A^.o, which is de- 
termined by the general formula (A12|. In the limit of 
As ^ |Aq| we can use the asymptotics of the function 
U{a) in the limit a ^ 1, U{a) « ln(4/a) — 1 leading to 
the following simple result 



A 



s,0 



lBa\^a\ 



In 



4|A^| 
A, 



1 



(A13) 



The sign of A^ g determines net effect of the s± super- 
conductor on the s superconductor (positive vs negative 
proximity). As one can expect, the gaps aligned with 
As enhance superconductivity and the gaps anti-aligned 
with As suppress superconductivity in the s supercon- 
ductor. The relative contributions are mostly determined 
by the electrical coupling between s-superconductor and 
the s± bands. 

Weak spatial dependence of As is determined by the 
components Ag.m with to > 0. At T — these compo- 
nents can be presented as 



As.m 



1 



Wra ds£,s 



E 



As- A^ 



(AM) 



with 



dz 



lo (z2 + 1)3/2 



1 + 



v/i^TA2 + /32 



(z2 + 1) ^z2 + (A„/As)2 + /3„ 



and Pm — [T^mS^s.^l ds) ■ The first integral can be actu- 
ally evaluated in general case 

1 



Wrr 



and it has asymptotics Wm ~ lii2/3„i for ^ 1. In the 
limits |Aq| As and Ag^m ^ |Aq,| the second integral 
can be evaluated as 



J Q 



dz 



1 



1 



Prn |A„| 

Collecting terms, we obtain 



v/^2 + (A^/A^)2z + /3„ 

1 , 2As/3™ 
m ■ 



As 



E 



1 



ln(As/|A„|) 



In 



2 {TTm^s.A/d-sY 



2ds/Cs A. - A^ 

(tTto)^ iBa 

(A15) 

Using relation (|a;|-l)2 = 1+4^^^^ ''°(iZ'f ^ ' can 
approximately present the gap correction in real space as 

■(.T + ds)' 1" 



As(.).As,-|E^^ 

?s „ IBa 

^ ln(As/|A„|) 



2d2 6 



2 In {TrV2^s.A/ds 



(A16) 



Correspondingly, for the Green's function for the same 
limits we derive 



1 As - A„ 



^ ^^2 + A2 jBc 



1 f X + da 



2 V 



(A17) 
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Typically, the spatial dependence of the order param- 
eter is correlated with the sign of proximity: As(a;) 
increases towards the boundary fo r pos itive proximity 
and vice versa. Analyzing Eq. (A16), however, we 
can conclude that this is not always the case. In- 
deed, the quadratic term changes sign roughly when 
(Ai - A,)/7bi - (A, + |A2|)7b2 = 0. At this point 
A,,o oc ln(CAi/A,)/Ai - In (CjAal/A,) /IA2I and it 
is negative for Ai > IA2I meaning that in some nar- 
row range of parameters As (a;) will increase towards the 
boundary in spite of negative proximity effect. 



The component ^a.A has to be found from the following 
equations 



^IuKa - *"^A = -A„ 



(A25) 
(A26) 



We can again find ^a,A{x) and Aa{x) using Fourier 
transform 



, COS qmX, 



A„ = ^ A„ 



m cos QniX 



2. s±-wave superconductor 

We can evaluate corrections to the s± gap parameter 
and Green's functions following the same general route. 
The first-order correction to with respect to the cou- 
pling strength 7^^ is determined by the following equa- 
tion and boundary conditions. 



— G $" 
Gs 



CaG^^a = —(A, 
7Sa 



Aq), at X 



(A18) 
(A19) 



with Ga ~ uj/\/uJ^ + A2 and = at x = -d±. The 
self-consistency condition for the linear correction 



2nT 



E 

P, ui>0 



(A20) 



can be rewritten as 



2'kT 



u>>Q 



A, 



$1 



InlA, 



E 



(A21) 



where A is the largest eigenvalue of the matrix X^p. 

Similar to the s-wave case, we can split $0 into the 
contributions induced by the boundary condition and by 
the correction to the gap parameter, = <ba.b + ^a,A- 
The equation and the boundary condition for ^a,b si's 



^2 (j)" A 



0, 



eaG„$U = - — (A.-A„) 



(A22) 
(A23) 



with = /?a/(2yZJ^+A|) = C^A\^^a\/V^o^+M 
and A = Da/{2\Aa\). The solution for $0 b{x) is given 
by 



G, 



(A, 



^ cosh[(a;-d±)/^nA 



with q„i = mi: /d±. From Eq. ( A25 1 we immediately find 

A 

(A27) 



1 £2 „2 



Substituting ^a,A hrto the gap equation, we obtain re- 
lation connecting A^ „ with '^p^b.mi Fourier components 
of ^p,b(x), 



0.b,m 



27tT J2 Um 



3/2 ' 



(A28) 



5^m,a — 2TrT ^ ^ 

CJ>0 . 



In- 



Fourier transformation of the result in Eq. ( A24 ) gives 

(2 - S,n) ec.,./ (rf±£a) + As - A„ 



s.b.m 



1 + £2 fl2 



(A29) 

It is convenient to introduce the degenerate matrix, 
Wai3 = ~ ^~^^ai3, W11W22 - W12W21 = 0. Explicitly, 
the elements of this matrix are given by 



Wii 



W22 



XI/4 + X12X21-X-/2 
dctA ' 



W12 



M2 

det A 



^A^/4 + Ai2A2i + A_/2 



det A 



with A_ = All ~ A22 and det A = A11A22 — A12A21. The 
matrix Um,ai3 — [wap — S„i^q(5q,^] ^ can be presented as 

1 



Urn. 



Dj 



W22-Sm,2 -W12 
-W21 Wii-Y.m.1 



(A30) 



Du — — Sm_2Wii — S„i,l'f^22 + Sm^iEm_2- 



The equations (A27l, (A28l, and (A29) provide full for- 



(A24) mal solution of the problem. 
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For convenient comparison with numerical calcula- 
tions, we present the above results in the reduced form 



Q:,A,7n 



2t U, 



7s o 



A2 v/oT^TAf 



This corresponds to equations in the reduced units pre- 
sented in the text. 



For the TRSB state the functions are com- 

plex and we select the following parametrization <^g = 
LuieaiOse^f' and — iut tan ^e^'^ " . Taking derivatives 



^ S 



cos2 9 a 



-f i(p'„tan6'„ e"^°, 



Appendix B: Boundary conditions in 
6*- parametrization 

In this appendix we derive the boundary conditions 
(BCs) in 0-parametrization, appropriate for the differ- 
ent SC states we consider. We start with the multiband 
generalization of the Kupriyanov-Lukichev BCs: 



„ la 



iBc 



-G,($,-$„). 



(Bla) 
(Bib) 



la 



C«^a = sin(0„-0,). 



(B2) 
(B3) 



substituting them in Eqs. (Bl) and separating real and 



imaginary parts, we obtain the boundary conditions for 
general complex case 



In the aligned case all $'s can be chosen real. Thus we 
have: 

$ct = tan 6*0, ^s—ijjt&n 6s, Ga—cos9a, Gg— cos 9s, 

where 6''s are also real. With those the BC at the inter- 
face for the aligned state simplify to: 



CaVa sin 0a = siu ds COS {(fis - fa)' 

IB a 



(B4a) 
(B4b) 



C^s = ^y^ [cos 6a sin 0s -I- COS sin6'Q sin (ipa — 'Ps)] , 

(B4c) 

Cs'^'s siTL6s = E -J— sin6a cos {ipa - fs) ■ (B4d) 
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